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ABSTRACT 

SN 2007bi is an extremely luminous Type Ic supernova. This supernova is thought to be 
evolved from a very massive star, and two possibilities have been proposed for the explosion 
mechanism. One possibility is a pair-instability supernova with an Afco ^ 100 M© CO core 
progenitor. Another possibility is a core-collapse supernova with Mco ^ 40 M©. We investi- 
gate the evolution of very massive stars with main-sequence mass Mms = 100 — 500 M© and 
Zq = 0.004, which is in the metallicity range of the host galaxy of SN 2007bi, to constrain the 
progenitor of SN 2007bi. The supernova type relating to the surface He abundance is also dis- 
cussed. The main-sequence mass of the progenitor exploding as a pair-instability supernova 
could be Mms ~ 515 — 575 M©. The minimum main-sequence mass could be 310 M© when 
uncertainties in the mass-loss rate are considered. A star with Mms 110 — 280 M© evolves 
to a CO star, appropriate for the core-collapse supernova of SN 2007bi. Arguments based on 
the probability of pair-instability and core-collapse supernovae favour the hypothesis that SN 
2007bi originated from a core-collapse supernova event. 

Key words: stars: evolution — stars: massive — stars: mass-loss — supernovae: individual: 
SN 2007bi — stars: Wolf-Rayet. 



1 INTRODUCTION 

SN 2007bi was found as an extremely luminous Type Ic supernova 
(SN Ic). Spectral ana lyses deduced the pro duction of 3.6 — 7.4 M0 
of radioactive ^''Ni (Gal - Yam et al.ll2009h . The metallicity of the 
host galaxy of SN 2007bi, which is a subluminous dwarf galaxy, 
has been observed to be 12 + l og(0/H) = 8.15 ± 0.15, corre- 
sponding to Z = 0.2 - 0.4 Zq jYoung et alj|2010h . This metallic- 
ity is similar to those of low-metallicity galaxies undergoing 7-ray 
bursts (GRBs) associated with SNe. 

This SN was proposed as a pair-instability (PI) SN from ob- 
servations of the light curve and spectral analyses. The final mass 
was estimated to be ~ 100 M©. On the other hand, the light curve 
was also fitted by th e energetic core-collapse (CC) explosion model 
jMoriva et al.ll20ld) . The estimated progenitor was a 43— M0 CO 
core. 40 Mq of ejecta, containing 6.1 M© of ^^Ni, were thought to 
be ejected. Although there is a difference in the rise period of the 
light curves of these models, there were no observations during the 
period. Hence, the explosion mechanism of SN 2007bi has not yet 
been clarified. Observational features such as the total ^''Ni mass 
and the SN type will provide constraints for the progenitor. 

The evolution of very massive stars connects presupernova 
structures and stars at the main-sequence (MS) stage. The final stel- 
lar mass, the mass of the CO core and the surface composition de- 
pend on MS mass and metallicity through burning processes and 
mass loss. Calculations of the evolution of very massive stars will 
provide the relationship of these features to MS mass. 
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The relationship between the amount of Ni produced in a 
SN and the CO core mass has been evaluated for PI and CC SN 
models. The mass range of the CO cores for PI SN progenitors 
has been evaluated as Mco ~ 60 — 130 M© dHeger & WooslevI 
I2OO2I : lUmeda & Nomotd 



!t3^ 

were also suggested bv IWaldmai] j2008l) . The amo unt of ^^Ni 
produ ced in a PI SN increases with MS mass (Hegei '^ WooslevI 
l2002l : IUmeda & Nomotdl2002h . The amount of ^'^Ni ejected from 
a CC SN is related to the CO core mass and the explosion energy 
dUmeda & Nomotdl2008h . A SN with an Afco ^ 35 Mq progeni- 
tor can produce more than 3 Mq of ^''Ni when the explosion energy 
is _Ekin = 3 X 10^^ erg. These relations connect the amount of ^''Ni 
deduced from the observations of SN 2007bi with the theoretical 
MS mass of a very massive star. 



2OO2I). Smaller progenitors for PI SNe 



In this Letter, we evaluate the mass range of very massive stars 
with initial metallicity of Zo — 0.004 (= 0.2 Zq) appropriate for 
the progenitor of SN 2007bi. We calculate the evolution of very 
massive stars with MS mass range Mms ~ 100 — 500 Mq taking 
into account uncertainties in mass-loss rate. We show the CO core 
mass as a function of the MS mass. Then, we constrain the MS and 
CO core masses using the amount of ^''Ni observed in SN 2007bi. 
We also discuss the relationship between the SN type and surface 
He abundance. Finally, we deduce the ranges of MS mass in PI and 
CC SN models appropriate for SN 2007bi and discuss the possible 
explosion mechanism. 
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2 PROGENITOR MODEL 
2.1 Stellar evolution model 

In order to estimate the progenitor mass of SN 2007bi, we cal- 
culated the evolution of very massive stars from H burning until 
central He exhaustion. We considered a mass range of zero-age 
MS stars of Mms = 100 - 500 M©. The initial metallicity was 
set to be Zq = 0.004 (— 0.2 Zq). The initial mass fractions of 
and ^He were set to be 0.7492 and 0.2468, respectively. The 
relative abundances of spe cies heavier than He were derived from 

I Anders & Grevessd j 19891) . 

The stellar evolution cod e was updated from lSaio et al.l i ll988h 
and lUmeda & Nomotd i2005l) . We used a nuclear reaction network 
consisting of 282 species of nuclei from n, p, to Br to calculate 
the time evolution of the chemical composition distribution and nu- 
clear energy generation. We adopted reaction rates from the current 
version of REACLIB ( Rauscher & Thielemann 200 0). The reactio n 
rate of ^^C(a, 7)^®0 was taken from Caughl an & Fowled jl988l) . 
multiplied by a factor of 1.6. For the convection treatment we as- 
sumed the Schwarzschild criterion for convective instability, and 
a diffusive treatment for the convective mixing. The opacity was 
evaluated using tab les of OPAL opacity ( Iglesias & Roge rs 1996), 
molecular opacity ^Ferguson et"aill2005h , and conductive opacity 
jCassisi et al]|2007h . 



2.2 Mass-loss rate 

Since mass-loss rates suffer from various uncertainties, we present 
models computed with the following three different mass-loss 
recipes. 



2.2.1 Case A (standard case) 

The mass-loss rate of MS stars with effective temperat ure Teff > 
12 OOP K and H mass fraction X ^ 0.4 was taken from lVink et al.l 
( l200lh . The mass-loss rate of stars wi th low surface temper ature 
Tcff s; 12 000 K was adopted from ' de Jager et all h988t) . We 
multiplied the mass-loss rate by the metallicity dependence fac- 
tor (.Z/Z(7))°'^"^. The power index 0.64 is the same value as that of 
B supergiants in Vink et al] bOOlh and is consistent with observa- 
tional scaling jMauron & Josseliij|201 ih . 

We set the mass-loss rate of Wolf-Rayet (WR) stars with 
Toff > 12 000 K and X < 0.4 as 



a = 0.60 and 0.40 for WN and WC stars, respectively. The metal- 
licity dependence of the mass-loss rate of WR stars in the Large 
Magellanic Cloud and Small Magellanic Cloud was observed as 

n0.8±0.2 WW ^ fV.. /r,\0.6±0.2 



(Zo/Z 



for WN stars and oc (Zo/Zq)" 



log(-M) = -11.0 + 1.29 log(L/L0) 

+ 1.731ogmin(Fs,0.98) +0.471ogmax(Z,0.02) 
+ alog(Zo/0.02) [log(M0yr-^)], (1) 

where L and Ys are the luminosity and helium mass fraction at 
the surface. The factor for the initial metallicity dependence a was 
set to be 0.86 and 0.66 for WN and WC stars, respectively. WR 
stars were classified into WN, WC and WO stars in accordance 
withlGeorgy et alj (l2009f) . We took into account equation (22) in 
'Nugis & Lamers ( 2000') and the initial-metallicity dependence in 
Vink & de Koter (^005). Extreme, episodic mass loss in WR s tars, 
which has been observed in some SNe (e.g. lFolev et"^l200l for 
SN 2006jc) was not taken into account in the present study. 



2.2.2 Case B (larger mass loss during the WR phase) 

We adopted the same mass-loss rates in the MS and low surface 
temperature phases, and rates with metallicity-dependent indices 



forWC 

stars jCrowthej2607l) . If the scaling index is small, the correspond- 
ing mass-loss rate is large. 

2.2.3 Case C (smaller mass loss during whole lifetime) 

We reduced the mass-loss rate in case A by a factor of 2 dur- 
ing the whole lifetime. Overestimations of the empirical mass-loss 
rate in OB stars and WR stars have bee n discussed (see the re- 
view in lCrowthe3l2007l : |Pulse et al.ll2008h . The reduction of mass- 
loss rate depends on wind clumping (e.g. iFuUerton et alj l2006l : 
Mokiem et al. 2007) and the degree of clumping has large uncer- 
tainties jPulse et alj|2008l) . It is considered that a reducti on wit hin 
a factor of 2 from the theoretical evalu ation cVink et al.ll200l | ) is 
allow ed from stellar evolution models jHirschi 20081 ; Pulse et afl 
[2008^. The effect of wind clumping in WR stars was taken into 
account by Nugis & Lamers ( 2000). The reduction of global mass- 
loss rate of WR stars by wind cl umping is a factor of 2 — 4 relative 
to the homogeneous model tCrowthej|2007h . The surface tempera- 
ture of very massive stars does not decrease to log Teff < 3.6, and 
dust-driven winds might be ineffective. In this case, wind clumping 
might be effective for low surface temperature. 



3 RESULTS 

3.1 Final mass of very massive stars 

Fig. 1(a) shows the final mass Mi and type of WR star as a function 
of the MS mass. The final progenitor mass depends on the mass- 
loss rate as well as on the MS mass. The dependence on the mass- 
loss rate is larger for more massive stellar models. In case A the 
final mass increases with the MS mass. These stars evolve to WO 
stars except for the one with Mms = 100 Mq. The final mass 
of the stars in case B rises with the MS mass except for the range 
of 170 ^ Mms < 250 M©. All of the stars evolve to WO stars. 
The increasing final mass is also shown in case C. The final type of 
the star depends on the MS mass. The relatively 'low'-mass stars 
with 100 ^ Mms ^ 120 Mq evolve to WN stars. The stars with 
130 ^ Mms ^ 140 Mq and 200 Mms ^ 300 Mq evolve to WC 
stars. The former WC stars lose their remaining He envelope during 
carbon burning. The He envelope of the latter stars has mostly been 
stripped away before carbon burning. The other stars evolve to WO 
stars. 

The mass of the CO core of an evolved star is important for 
the explosion mechanism of SN 2007bi. Fig. 1(b) shows the CO 
core mass Afco as a function of the MS mass. The mass of the CO 
core is defined as the large st mass coordinate satisfy ing a He mass 
fraction smaller than 10"^ jUmeda & Nomotol200a) . The CO core 
mass of a star with a given Mms and mass-loss rate is smaller by 
~ 10 — 13 per cent than the corresponding final mass. 



3.2 Ejected amount of ^''Ni and mass of CO core 

From analyses of the light curve and spectra of SN 2007bi, the ejec- 
tion of 3.6 — 7.4 Mq of ^^Ni has been derived. The average amount 
of ^''Ni is 5.3 M0. The amount of ejected ^''Ni has been connected 
to the CO core mass in previous studies. Here we constrain the 
range of MS mass from the amount of ^^Ni in SN 2007bi. 
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Figure 1. The final mass (a) and CO core mass (b) as a function of the 
MS mass. Solid, dashed, and dotted lines correspond to cases A, B and C, 
respectively. Squares, triangles, crosses and circles indicate WN, 'He-rich' 
WC (see Section 3.3), WC and WO stars. Dark and light shaded regions de- 
note the mass ranges of the CO cores appropriate for SN 2007bi to explode 
as a PI SN and a CC SN, respectively. 



Nucleosynthesis studies of PI SNe have indicated that PI SN 
models with 95 < A/qq < 105 Mq produce 3 - 10 Mr 



Heger & Wooslevl2002l) . The PI SN models of lUmeda & Nomotd 
2002h obtained a similar result. Therefore the progenitor of the 



PI SN corresponding to SN 2007bi should be a CO core with 95 
< Mco < 105 Mq. 

The MS mass range deduced from the CO core mass is shown 
as the dark shaded region in Fig. 1(b). This criterion is not satisfied 
in the MS mass range of Mms ^ 50QMq in cases A and B. The 
MS mass range in case A is estimated to be 515 ^ Mms ^ 575 Mq 
from linear extrapolation in Fig. 1(6). In case C, the appropriate 
range reduces to 310 ^ Mms ^ 350 Mq. 

The amount of ^®Ni produced in a CC SN model gives a 
lo wer limit on the progenit or mass. We expect from the result 
of lUmeda&Nomotol ( l2008h that a SN with a kinetic energy of 
~ 3 X 10 erg can produce more than 3 M© of ^'^Ni if the CO 
core of the progenitor is larger than ~ 35 Mq. The upper limit 
of the MS mass might be the lowest mass of a PI SN progeni- 
tor. Theoretical studies of PI SNe have indicated that a CO core 
larger than > 60 Mq exp lodes as a PI SN (Heger & Wooslev 2002; 
lUmeda & Nomotoll20o3) . We consider the upper limit of a CC SN 
progenitor as ~ 60 Mq . Therefore, the CO core mass appropriate 
for explaining SN 2007bi with a CC explosion is 35 < Mco ^ 60 

Mq. 



The range of MS mass for CC SN models is shown as the light 
shaded region in Fig. 1(b). The range extends from 100 to 280 Mq 
in case A. In case B, all models except for Mms = 100 Mq will 
explode as CC SNe appropriate for SN 2007bi. On the other hand, 
in case C, the mass range is limited to 100 < Alms ^ 170 Mq. 



3.3 Surface He abundance for SNe Ib/Ic 

SNe Ic are characterized by weak or absent He spectra. However, 
a quantitative criterion to distinguish between SNe lb and Ic has 
not been theoretically established. The He lines are considered to 
appear because of the excitation of He by non-ther mal electron s 
excited by 7-rays from the decays of ^^Ni and ^^Co ( lLucvl[l99lh . 
The strength of the He lines should be sensitive to the amounts of 
He and ''"Ni, the amount of matter in intermediate layers between 
Ni and He layers which attenuate the 7-rays, the degree of mix- 
ing of Ni into He layer, etc. Criteria to distinguish between SNe 
lb and Ic were discussed using the total He m ass of a progenitor 
dWellstein & Langeil [l999l : lOeorgy et alj|2009l : Yoon et al. 20 13) 
or the He mass fraction at the outermost layers (Yoon et al. 2010^ 
The effects of thickness of the int ermediate layers and the d egree 
of mixing have been investigated jWooslev & Eastmanlll997h . We 
discuss the possibility of a SN Ic progenitor by considering the to- 
tal He mass, the He mass fraction at the surface and the mass ratio 
of He to the intermediate layers. 

The first criterion is based on the total He mass. We consider 
two cases for the He mass limit of SNe Ic: 0.5 and 1.5 Mq. In 
previous studies, the He mass limit o f SNe Ic w as assumed to be 
0.5 MQjWellstein&Langer 1999; Y oon et alJb olO) or 0^6 Mq 
dOeorgy et alj[2009l) . On the other hand, Georg y et alj ^200^ re- 
ported that the choice of the total He mass limit betw een 0.6 and 
1.5 M q hardly affects the Mms ranges for SNe Ib/Ic. lYoon et all 
( l2010l) suggested in discussion that He lines are not seen in early- 
time spectra even though the total He mass is as large as 1.0 Mq 
if He is well mixed with CO material having Ys < 0.5. Fig. 2(a) 
shows the total He mass versus the MS mass. When the He mass 
limit is 0.5 Mq, the MS mass is limited in the very narrow range 
of 110 — 120 Mq in case A and 100 — 115 Mq in case B. In case 
C no progenitors explode as SNe Ic. When the He mass limit is 
1.5 Mq, all progenitors except for WN stars, an 'He-rich' WC star 
with Mms = 130 Mq in case C, and stars with Mms > 350 Mq 
in case C will explode as SNe Ic. 

The second criterion is the He mass fraction at the surface. We 
set this criterion a s Ya ^ 0.5 in accordance with the discussion in 
lYoon et alj ilOld) as mentioned above. Fig. 2(b) presents the He 
mass fraction at the surface as a function of the MS mass. All WO 
stars and WC stars except for the He-rich WC will explode as SNe 
Ic. 

We also consider the mass ratio of He to the intermediate lay- 
ers. The mass of the intermediate layers is assumed to be Mint = 
M{ — 9.0 Mq. The mass of 9.0 Mq corresponds to the amount 
of ejected ^''Ni and the central remnant in the CC SN model. The 
masses of ^^Ni and the central remnant were evalu ated to be 6.1 
and 3 .0 Mq, respectively, for the 43— Mq progenitor jMoriva et al.l 
I2OIOI) . Although the amount of '^'^Ni would be smaller than 9 Mq 
and there is no remnant in the PI SN model, the difference is not 
important because the final mass of the PI SN progenitors is ~ 100 
Mq. The mass ratio of He to the intermediate layers is shown in 
Fig. 2(c). The ratios of WN stars and the He-rich WC star are 
clearly larger than those of WO and the other WC stars. We may set 
the criterion of the ratio as ^ 0.02 — 0.03. In this case, the range of 
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Table 1. The range of MS mass for PI and CC SN models appropriate for 
SN 2007bi deduced from the mass range of the CO core and the surface He 
abundance. The ratio of the probability of explosion as a PI SN to explosion 
as a CC SN for SN 2007bi, rpj/cc. is also listed (see Section 4). 
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Figure 2. The total He mass in the envelope (a), the He mass fraction at the 
surface (b), and the mass ratio of He to the intermediate layers (c) plotted 
against the MS mass. SoHd, dashed and dotted Hnes correspond to cases A, 
B and C, respectively. Squares, triangles, crosses and circles indicate WN, 
'He-rich' WC, WC and WO stars. In panel (a), dark and light shaded regions 
satisfy the criteria A/(Hc) ^ 0.5 and 1.5 Mq, respectively. In panel (b), 
the shaded region satisfies the criterion Yb ^5 0.5. 



range of CO core mass and the surface He abundance. Table 1 
presents the MS mass range with different conditions of He abun- 
dance and mass-loss rate. The MS mass range strongly depends on 
the surface He conditions. If He lines appear with a total He mass 
larger than 0.5 Mq, then there is no mass range of the PI SN model 
for SN 2007bi. The MS mass range of the CC SN model is also 
limited to ~ 110 — 120 Mq. However , this condition may be too 
strict, as discussed bv lYoonetal1 ( l2010l) . If the progenitors explode 
as Type Ic SNe with M(He) ^ 1.5 Mq or ^ 0.5, the MS mass 
range is mainly determined by the condition of the CO core mass. 
Thus the minimum MS mass for the PI SN model is 310 Mq. The 
MS mass range of the CC SN model in case A becomes slightly nar- 
rower: 1 10 ^ Mms ^ 280 Mq. In cases B and C the mass range is 
determined by the mass of the CO core. 

We now discuss the probabilities of SN 2007bi exploding as a 
PI SN and a CC SN, considering the MS mass range and Salpeter 
initial mass function (IMF) (dN/dMus oc M^l'^'^). The ratio 
of the probability of explosion as a PI SN to explosion as a CC 
SN, rpi/cC' is also listed in Table 1. If He lines are not seen for 
Af(He) s; 1.5 Mq or Ys ^ 0.5, the probability of the CC SN 
model is about 40 times larger than that of the PI SN model in 
case A. In case B, the PI SN model is not adequate for SN 2007bi 
because of large mass loss. In case C, on the other hand, the prob- 
ability of the CC SN model is larger than that of the PI SN model 
by a factor of 5. However, if He lines are observed with a total He 
mass of A/ (He) > 0.5 Mq, the PI SN model is inadequate for 
SN 2007bi in all mass-loss cases. It is important to establish more 
definite classification criteria for SNe Ic. 

It has been argued that very massive CO cores would collapse 
to for m black holes and become dark or very faint SNe (e.g., re- 
view m ISmarttl2009b . In this case, the probability of explosion as a 
CC SN would become much smaller. Observational or theoretical 
evaluation of the fraction of dark or faint SNe in the death of very 
massive stars will more precisely constrain the probability of the 
CC SN model. 



MS mass is the same as that deduced from the criterion of the He 
mass fraction at the surface. 

4 RANGE OF MAIN-SEQUENCE MASS FOR THE SN 
2007BI PROGENITOR 

We evaluate the ranges of MS mass appropriate for explaining SN 
2007bi by the PI SN and CC SN models taking into account the 



5 CONCLUDING REMARKS 

We have investigated the evolution of very massive stars with 
Mms = 100 - 500 Mq and Zo = 0.004 from H burning to central 
He exhaustion to constrain the progenitor of an extremely luminous 
Type Ic SN 2007bi. If SN 2007bi is a PI SN, the progenitor could 
be evolved from a MS star with 515 < Mms < 575 Mq. If SN 
2007bi is a CC SN, the appropriate range of the MS mass is 1 10 < 
Mms ^ 280 Mq. When we take into account a Salpeter IMF, the 



A progenitor for SN 2007bi 5 



probability of SN 2007bi exploding as a CC SN is about 40 times 
larger than for a PI SN. If the mass-loss rate is small, the minimum 
mass for the progenitor becoming a PI SN is about 310 M©. The 
mass range changes to Mms ~ 310 — 350 Mq for the PI SN model 
and Mms 135 - 170 Mq for the CC SN model. If the mass-loss 
rate is large or He lines appear for an He mass of M(He) > 0.5 
Mo, then SN 2007bi should have exploded as a CC SN. Although 
the light curve of SN 2007bi was fitted by both the PI SN and CC 
SN models, the CC SN model is favoured as explaining the explo- 
sion of SN 2007bi from arguments based on the probability ratio of 
the PI SN model appropriate for SN 2007bi to the CC SN model. 

We should note that the conditions of He abundance for the 
progenitor to become a Type Ic SN strongly affect the probability 
of a PI SN. If a total He mass M(He) > 0.5 Mq enables the 
formation of He spectra, the PI SN model would be inadequate for 
SN 2007bi even for a small mass-loss rate. It is quite important to 
evaluate more definite criteria to classify SNe Ic. We also note the 
probability of direct collapse without bright SNe for CC SN cases. 
This would reduce the probability of a CC SN for SN 2007bi. It is 
necessary to evaluate theoretically or observationally the fraction 
of black hole formation without bright SNe in CC SNe. 

If the progenitor of SN 2007bi or a more massive star was ro- 
tating very fast, it could explode as a SN Ic associated with a GRB. 
We may not have seen such a GRB-SN because of the orientation 
effect. This is also interesting point to investigate in the future. 



Lucy, L. B., 1991, ApJ, 383, 308 

Mauron, N., Josselin, E., 2011, A&A, 526, A156 

Moriya, T, Tominaga, N., Tanaka, M., Maeda, K., Nomoto, K., 

2010, ApJL, 717, L83 
Mokiem, M. R. et al., 2007, A&A, 473, 603 
Nugis, T, Lamers, H. J. G. L. M., 2000, A&A, 360, 227 
Pulse, J., Vink, J. S., Najarro, R, 2008, A&AR, 16, 209 
Rauscher, T, Thielemann, F.-K., 2000, At. Data Nucl. Data Ta- 
bles, 75, 1 

Saio, H., Nomoto, K., Kato, M., 1988, ApJ, 331, 388 
Smartt, S. J., 2009, ARA&A, 47, 63 
Umeda, H., Nomoto, K., 2002, ApJ, 565, 385 
Umeda, H., Nomoto, K., 2005, ApJ, 619, 427 
Umeda, H., Nomoto, K., 2008, ApJ, 678, 1014 
Vink, J. S., de Koter, A., 2005, A&A, 442, 587 
Vink, J. S., de Koter, A., Lamers, H. J. G. L. M., 2001, A&A, 369, 
574 

Waldman, R., 2008, ApJ, 685, 1103 

Wellstein, S., Langer, N., 1999, A&A, 350, 148 

Woosley, S. E., Eastman, R. G., 1997, In: Thermonuclear Super- 
novae, eds. Ruiz-LaPuente P, Canal, R., & Isem, J., NATO ASI 
Ser. C, Vol. 486, Kluwer Publishers, 821 

Yoon, S.-C, Woosley, S. E., Langer, N., 2010, ApJ, 725, 940 

Young, D. R. et al., 2010, A&A, 512, A70 



ACKNOWLEDGMENTS 

We acknowledge the anonymous referee for many valuable com- 
ments. We thank Hideyuki Saio for providing the stellar evolution 
code and for useful comments. We are grateful to Masaomi Tanaka, 
Nobuyuki Iwamoto, and Ken'ichi Nomoto for useful discussions. 
This work was supported by the grants-in-aid for Scientific Re- 
search (20041005, 20105004) from the MEXT of Japan. 



REFERENCES 

Anders, E., Grevesse, N., 1989, Geochim. Cosmochim. Acta, 53, 
197 

Cassisi, S., Potekhin, A. Y., Pietrinfemi, A., Catelan, M., Salaris, 

M., 2007, ApJ, 661, 1094 
Caughlan, G. A., & Fowler, W. A., 1988, At. Data Nucl. Data 

Tables, 40, 283 
Crowther, R A., 2007, ARA&A, 45, 177 

de Jager, C, Nieuwenhuijen, H., van der Hucht, K. A., 1988, 
ApJS, 369, 574 

Ferguson, J. W, Alexander, D. R., Allard, R, Barman, T, Bod- 
narik, J. G., Hauschildt, P. H., Heffner-Wong, A., Tamanai, A., 
2005, ApJ, 623, 585 

Foley, R. J., Smith, N., Ganeshalingam, M., Li, W., Chomock, R., 
Filippenko, A. V., 2007, ApJL, 657, L105 

Fullerton, A. W, Massa, D. L., Prinja, R. K., 2006, ApJ, 637, 1025 

Gal- Yam, A. et al., 2009, Nat, 462, 624 

Georgy, C, Meynet, G., Walder, R., FoUni, D., Maeder, A., 2009, 
A&A, 502,611 

Heger, A., Woosley, S. E., 2002, ApJ, 567, 532 

Hirschi, R., 2008, in Clumping in Hot Star Winds, eds. W.-R. 
Hamann, A. Feldmeier, L. M. Oskinova, Universitatsverlag Pots- 
dam, Potsdam, p. 9 

Iglesias, C. A., Rogers, R J., 1996, ApJ, 464, 943 



